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ABSTRACT 

We investigate the alignment of the spin of dark matter halos relative (i) to the sur- 
rounding large-scale filamentary structure, and (ii) to the tidal tensor eigenvectors 
using the Horizon 4/1 dark matter simulation which resolves over 43 million dark 
matter halos at redshift zero. We detect a clear mass transition: the spin of dark mat- 
ter halos above a critical mass Mq w 5(±1) ■ tends to be perpendicular to 
the closest large scale filament (with an excess probability up to 12%), and aligned 
with the intermediate axis of the tidal tensor (with an excess probability of up to 
40%), whereas the spin of low-mass halos is more likely to be aligned with the clos- 
est filament (with an excess probability up to 15%). Furthermore, this critical mass 
is redshift-dependent, scaling as M^^.jj(z) w Mg ■ (1 + z)"'>'-' with 7s = 2.5 ± 0.2. 
A similar fit for the redshift evolution of the tidal tensor transition mass yields 
Mq ~ 8(±2) ■ IO^^Mq and 7f = 3 ± 0.3. This critical mass also varies weakly with 
the scale defining filaments. 

We propose an interpretation of this signal in terms of large-scale cosmic flows. 
In this picture, most low-mass halos are formed through the winding of flows embed- 
ded in misaligned walls; hence they acquire a spin parallel to the axis of the resulting 
filaments forming at the intersection of these walls. On the other hand, more massive 
halos are typically the products of later mergers along such filaments, and thus they 
acquire a spin perpendicular to this direction when their orbital angular momentum 
is converted into spin. We show that this scenario is consistent with both the mea- 
sured excess probabilities of alignment w.r.t. the eigen-directions of the tidal tensor, 
and halo merger histories. On a more qualitative level, it also seems compatible with 
3D visualization of the structure of the cosmic web as traced by "smoothed" dark 
matter simulations or gas tracer particles. Finally, it provides extra support to the 
disc forming paradigm presented by Pichon et al. (2011) as it extends it by character- 
izing the geometry of secondary infall at high redshift. 

Key words: large-scale structure - galaxies: halos - galaxies: formation - method: 
numerical. 



1 INTRODUCTION 

Over the past decades, numerical simulations and large red- 
shift surveys have highlighted the large-scale structure of 
our Universe (LSS), a cosmic web formed by voids, sheets, 
elongated filaments and clusters at their nodes (Bond et al. 
1996). This structure is believed to be the result of the lin- 
ear growth of primordial gaussian fluctuations in a nearly 
homogeneous density field followed by the non-linear col- 
lapse of overdense regions into dark matter halos which 
then accrete mass and merge as described by the hierarchi- 



cal model. The current paradigm of galaxy formation states 
that collapsing protogalaxies acquire their spin (ie their an- 
gular momentum) by tidal torquing because of a misalign- 
ment between their inertia tensor and the local gravitational 
tidal tensor at the time of maximum expansion; this is the 
basis of the so-called Tidal Torque Theory (TTT hereafter, 
Hoyle 1949; Peebles 1969; Doroshkevich 1970; White 1984; 
Catelan & Theims 1996; Crittenden et al. 2001; Schafer 2009, 
for a recent review). According to this theory, the spin direc- 
tion should initially be correlated with the principal axes of 
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the local tidal tensor, defined as the traceless part of the Hes- 
sian matrix of the gravitational potential field. One therefore 
expects to detect correlations between actual galactic angular 
momenta and large-scale structures if non-linear processes 
have not modified their direction. 

For many years, both observers and theorists have thus 
endeavoured to detect these correlations in real surveys and 
cosmological N-body simulations. Nonetheless, the results 
remain in part contradictory because of the lack in resolution 
together with the difficulty of properly defining large-scale 
filamentary structures. 

For instance, using N-body simulations Hahn et al. 
(2007b), Sousbie et al. (2009) and Zhang et al. (2009) found 
that halo spins are preferentially oriented perpendicularly to 
the fUaments whatever their mass, with Faltenbacher et al. 
(2002) measuring a random distribution of the spin in the 
plane perpendicular to the filaments, while Hatton & Ninin 
(2001) claimed an alignment between spin and filament. 
More recently a consensus seemed to have emerged when 
several works (Bailin & Steinmetz 2005; Aragon-Calvo et al. 
2007; Hahn et al. 2007a; Paz et al. 2008) reported that large- 
scale structure - filaments and sheets - influenced the di- 
rection of the angular momenta of dark matter halos in 
a way originally predicted by Sugerman et al. (2000) and 
Lee & Pen (2000). These studies pointed towards the detec- 
tion of a mass-dependence orientation of the spin, arguing 
for the first time that the spin of high-mass halos tends to 
lie perpendicular to their host filament, whereas low-mass 
halos have a spin preferentially aligned with it. Neverthe- 
less, the detected correlation remains weak and noisy and no 
full explanation for these findings were highlighted except 
e.g. Bailin & Steinmetz (2005) who suggested that the spin 
direction of the cluster and group mass halos (as opposed to 
galaxy mass halos) could come from mergers along the fila- 
ments. 

However, Hahn et al. (2010) repeated this measurement 
in a cosmological hydrodynamical re-simulation of one 
large-scale cosmic filament and found a different result: the 
spin of high-mass halos is aligned with the filament; the spin 
of low-mass halos is along the intermediate eigen-direction 
of the tidal tensor in low-density regions and along the third 
eigen-direction (i.e neither the intermediate direction nor the 
filament's) in higher density regions at higher redshift (z=l); 
finally no signal exists for high-density region at low redshift. 

Beyond numerical simulations, Lee & Erdogdu (2007); 
Trujillo et al. (2006); Navarro et al. (2004); Flin & Godlowski 
(1990, 1986); Godlowski & Flin (2010) foimd correlations in 
observations between the rotational axis of galaxies and the 
surrounding large-scale structures (e.g voids and local tidal 
shears) unlike Dekel (1985) who did not find any correlation. 
To be more specific, Flin & Godlowski (1986, 1990) first dis- 
covered that the spin of galaxies was not isotropically ori- 
ented w.r.t the Local Supercluster plane but more likely to 
be aligned with it; Navarro et al. (2004) confirmed this ob- 
servation; Trujillo et al. (2006) found that in the SDSS and 
2dFGRS, the rotational axis of the spiral galaxies located in 
the walls surrounding voids lie preferentially in the plane of 
these voids; Lee & Erdogdu (2007) analyzed the galaxies of 
the 2Mass Redshift Survey and found correlations between 
their spin and the local tidal tensor and Godlowski & Flin 
(2010) focused on the galaxy groups in the Local Superclus- 
ter and found correlations in their orientation suggesting that 



the two brightest galaxies and then the galaxy groups were 
hierarchically formed in the same filament with their ma- 
jor vector aligned with this host filament . The results of 
Navarro et al. (2004); Trujillo et al. (2006) support the predic- 
tions of TTT (Lee & Pen 2000) namely that, assuming that 
the inertia and tidal tensors are uncorrelated, galaxies's spin 
should be preferentially aligned with the intermediate eigen- 
direction of the tidal tensor (in particular in the plane of 
the voids). Recently, however, Slosar & White (2009) claimed 
that in contrast to previous studies, they found no departure 
from randomness in the SDSS while studying the orienta- 
tion of the galaxy spin with regards to the voids in which 
they are located. The methods used in the latter study has 
been improved by Varela et al. (2011) who used the SDSS 
(DR 7) and morphological classification from the Galaxy Zoo 
Project and found that galaxy disks are more likely to lie 
in the plane of voids i.e. their spin tend to be perpendic- 
ular to the void they are located in, which seems in dis- 
agreement with Lee & Erdogdu (2007); Trujillo et al. (2006); 
Navarro et al. (2004). 



In short, even if one can claim that a trend is slowly 
emerging, quantitative evidence of spin alignment with the 
filaments and tidal tensor eigen-directions remains at this 
stage weak and somewhat inconclusive. Hence, in this pa- 
per, we propose to revisit the issue and quantify the align- 
ment between the spin of dark matter halos and the fila- 
mentary pattern in which they are embedded (together with 
the alignment between the spin and the tidal tensor princi- 
pal axes) using a very efficient topological tool, the skele- 
ton (Sousbie et al. 2009). This tool provides a robust and 
mathematically well defined reconstruction of the cosmic 
web filaments. We apply it to the Horizon in simulation 
(Teyssier et al. 2009), a 2 h^^Gpc on a side cubic volume of 
the universe containing over 67 billion dark matter particles 
which provides an unprecedented catalog of 43 million dark 
matter halos with masses > 2 • 10^^M(3. We then interpret 
our results in the framework of the dynamics of large-scale 
cosmic flows. 



Section 2 briefly presents the Horizon in simulation and 
the topological tool implemented to identify the loci and 
orientation of filaments. It then reports the correlations de- 
tected between the orientation of the spin of dark matter ha- 
los and filaments and its redshift evolution. Section 3 is de- 
voted to the physical processes which induce these correla- 
tions. It also illustrates them using dark matter halo merging 
histories, smoothed dark matter simulations and hydrody- 
namical simulations. Section 4 provides conclusions and dis- 
cusses prospects for our understanding of galaxy formation 
within its cosmic environment. Appendix A gives the cor- 
relations measured between spin direction and tidal eigen- 
directions, which are in agreement with the cosmic dynamics 
arguments of Section 2. Appendix B presents a visual quanti- 
tative estimation of the spin of the circum-galactic medium. 
Appendix C sums up all the tests we have performed to as- 
sess the robustness of the correlations presented in Section 2. 
Finally, Appendix D presents the dependence of the transi- 
tion masses with the smoothing length and the non-linear 
mass as a function redshift. 
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2 SPIN-FILAMENT CORRELATIONS 

Let us first account for the robust correlation between the 
dark matter halo's spin and the orientation of the filaments 
of the cosmic web. 



2.1 Virtual data sets 

This study uses the Horizon in N-body simulation 
(Teyssier et al. 2009) which contains 4096^ dark matter parti- 
cles distributed in a 2 /i^^Gpc periodic box to investigate the 
spin alignment of dark matter halos relative to their large- 
scale structure environment. This simulation is characterized 
by the following ACDM cosmology: Qm = 0.24, Oa = 0.76, 
n = 0.958, Ho = 73 km-s^i-Mpc^^ and = 0.77 within one 
standard deviation of WMAP3 results (Spergel et al. 2003). 
These initial conditions were evolved non-linearly down 
to redshift zero using the AMR code RAMSES (Teyssier 
2002), on a 4096^ grid. The motion of the particles was fol- 
lowed with a multi grid Particle-Mesh Poisson solver using 
a Cloud-ln-Cell interpolation algorithm to assign these par- 
ticles to the grid (the refinement strategy of 40 particles as a 
threshold for refinement allowed us to reach a constant phys- 
ical resolution of 10 kpc, see the above mentioned two refer- 
ences). 

The Friend-of-Friend Algorithm (hereafter EOF, 
Huchra & Geller 1982) was used over 18^ overlapping sub- 
sets of the simulation with a linking length of 0.2 times the 
mean inter-particular distance to define dark matter halos. 
In the present work, we only consider halos with more than 
40 particles, which corresponds to a minimum halo mass 
of 3 ■ 10" M© (the particle mass is 7.7 ■ lO^M©). The mass 
dynamical range of this simulation spans about 5 decades. 
Overall, 43 million halos were detected at redshift zero (see 
fig. 1). This simulation was complemented by smaller (1024^ 
particles, boxsize 200 h^^Mpc and several 256^ particles, 
boxsize 50 /i^^Mpc leading to a particle mass of 6.2 ■ lO^M.j,) 
dark matter only simulations to address resolution issues 
(see Appendix C) and interpret the redshift evolution of the 
signal (see Section 2.3). 

Several topological techniques (Barrow et al. 1985; 
Stoica et al. 2005; Novikov et al. 2006) have been proposed to 
identify the complex cosmic network made of large voids, 
sheet-like structures and elongated filaments. These tech- 
niques rely on giving a mathematical definition (and a de- 
tection algorithm) of the filamentary pattern that our eye 
easily detects in the simulations. One recent criterion to 
classify structures as clusters, filaments, sheets or voids, 
uses the number of positive eigenvalues of the Hessian ma- 
trix of the density or potential fields (Hahn et al. 2007a; 
Aragon-Calvo et al. 2007; Forero-Romero et al. 2009). An- 
other interesting approach was followed by Platen et al. 
(2007, 2008) who used Watershed transforms to identify 
voids, filaments and walls. More recently Gonzalez &c Padilla 
(2010) introduced a method which relies on positions and 
masses of dark matter halos, and Bond et al. (2010) used Hes- 
sian eigen-directions to detect filaments. In this paper, we use 
the (global) 3D skeleton introduced by Sousbie et al. (2009). 
The underlying algorithm is based on Morse theory and de- 
fines the skeleton as the set of critical lines joining the max- 
ima of the density field through saddle points following the 
gradient. In practice Sousbie et al. (2009) define the peak and 




Figure 1. An 80 Mpc/h slice of the Horizon 4/7 simulation at red- 
shift zero. The box-size is 2Gpc/?! across. On top of the dark matter 
log-density (colour coded in levels of blue, from dark to light), all 
halos in that slice whose mass is larger than 3 ■ lO^^M:; are shown as 
yellow dots. As expected, these halos fall on top of the filamentary 
structure of the cosmic web. 




Figure 2. A 3D view of the skeleton of the Horizon 47r simulation 
measured from the dark matter distribution. The size of the box is 
2 Gpc/h. This paper analyses the relative orientation of the spin of 
the 43 million dark halos relative to this cosmic web. 
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void patches of the density field as the set of points converg- 
ing to a specific local maximum/ minimum while following 
the field lines in the direction/opposite direction of the gra- 
dient. The skeleton is then the set of intersection of the void 
patches i.e. the subset of critical lines connecting the saddle 
points and the local maxima of a density field and parallel to 
the gradient of the field. 

For this work, the ~70 billion particles of the Horizon in 
were sampled on a 2048^ cartesian grid and the density field 
was smoothed over 5 sigmas using mpsmooth (Prunet et al. 
2008), corresponding to a scale of 5 Mpc/h and a mass of 
1.9 X 10"M(.:,. Hence we are focusing on the large-scale struc- 
ture of the cosmic web. The corresponding cube was then 
divided into 6^ overlapping sub-cubes (with a buffer zone of 
100 voxels in each direction, large enough to cover the largest 
peak patches of the simulation, see Sousbie et al. 2009), and 
the skeleton was computed for each of these sub-cubes. It 
was then reconnected across the entire simulation volume to 
produce a catalog of segments which locally defines the di- 
rection of the skeleton. This skeleton is shown in fig. 2. Note 
that this skeleton (i.e what we will call filaments in the rest of 
this paper) depends on the choice we made for the smooth- 
ing length (5 Mpc//i). Appendix Dl investigates the effect of 
probing smaller smoothing scale on other sets of simulation. 

The hydrodynamical simulations used in this paper are 
described in Appendix B. 



2.2 Correlations between spin and filament axis 

In order to study the alignment between the spin of halos 
and the filamentary features of the cosmic web, we compute 
the skeleton of the LSS for the density field smoothed with 
the above quoted Gaussian scale R = 5 Mpc/fc which cor- 
responds to <j(-R) = 0.66. Thus we are considering the fil- 
aments that are mildly non-linear large-scale structures at 
cluster scales. 

In this paper, the spin of a given halo is defined as the 
following sum over its particles i: mpY^j{i[ — r) x (vj — v) 
where f is the center of mass of the FOF and v its mean ve- 
locity. We search for the five dark matter halos (regardless of 
their mass) closest to each filament segment (see Appendix C 
for alternative choices). We then measure the angle between 
the angular momentum of these halos and the direction of 
the filament segment and estimate the probability distribu- 
tion (PDF) of the absolute value of the cosine of this angle; 
this PDF, 1 + ^, measures the excess probability of alignment 
between the halo spin and the direction of the filament (note 
in particular that it is normalized for cos 9 between and 1; 
for aesthetic purpose only, data are symmetrically plotted for 
COS0 between —1 and 1; Appendix C briefly discusses the 
associated biases). The data is split by halo masses ranging 
from galactic to cluster masses and is displayed in fig. 3, the 
main result of this paper. 

A clear signal is detected. The orientation of the halo 
spin depends on the local anisotropy of the cosmic web, 
and on the dark matter halo mass: the spin of dark mat- 
ter halos is preferably perpendicular to their host filament 
at high mass (with an excess probability reaching 12 %), 
but turns into being aligned with the nearest filament di- 
rection at lower masses (with an excess probability of 15 %). 
This "phase-transition" is found to occur at M^^^^^{z = 0) ~ 




\ 



Figure 4. A very small subset of skeletons at different redshifts ex- 
tracted from the Horizon 4n simulation (see fig. 2), together with 
unit vectors showing the orientation of the spin of the correspond- 
ing dark matter halo with a mass above the transition mass. The spin 
is indeed perpendicular to the fUament for these massive halos. 

4(±1) • lO^^M;::, where M^j.;f is defined as the halo mass for 
which (cos 9) = 0.5. Fig. 4 shows an example segment of 
the large-scale filamentary network together with orientation 
of spins of massive haloes that graphically demonstrates for 
them the effect of spin-filament anti-alignment. Several san- 
ity checks have been carried out to assess the robustness of 
this signal and are summed up in Appendix C. 

These measurements of the spin-filament corre- 
lation trend confirm the previous results obtained by 
Bailin & Steinmetz (2005); Aragon-Calvo et al. (2007); 
Hahn et al. (2007a); Paz et al. (2008) with significantly im- 
proved statistics which allows us to quantify here the mass 
transition. 



2.3 Redshift dependence of transition mass 

The redshift-dependence of the transition mass was then in- 
vestigated on a set of smaller ACDM simulation (1024^ par- 
ticles in 200 h^^Mpc periodic boxes and 256^ particles in 50 
Jz^^Mpc periodic boxes). At high redshift we define the fil- 
amentary structure at the smoothing scale, R{z), chosen to 
maintain the same level of non-linearity as we had at redshift 
zero for Rg = 5Mpc//!. Thus, the smoothing scale R{z) is ob- 
tained from the implicit condition £7^ (_R(z),z) = a-^{RQ,Q), 
see Appendix D2. The halo spins continue to exhibit a transi- 
tion from alignment with the nearest filament at low mass to 
anti-alignment at high mass. The critical mass for the transi- 
tion, M^j.jj(z), is found to decrease with redshift as a power- 
law of z. Namely, 

M^rit « -^0(1 + z)"'^% 7s = 2.5 ± 0.2 , M'o~5- IO^^Mq. (1) 

This result is presented in Fig. 3 over the studied range 
z = — 4. Measuring the dependence of the z = transi- 
tion mass Mg on scale Rq (see Fig. Dl), we find some weak 
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Figure 3. Left; excess probability of alignment between the spin and the direction of the closest filament as measured from the 43 millions halos 
of the Horizon in simulation at redshift zero. Different colors correspond to different mass bins from 10^^(red) to lO^" M,3 (blue) as labeled. 
Thanks to the very large number of halos in each mass bin, the excess probability is quite well sampled and displays a clear departure from 
a uniform distribution. A transition mass is detected at Mjj = M^j.jj(z = 0) ~ 5(±1) ■ lO^^M,;;,: for halos with M > Mg, the spin is more 
likely to be perpendicular to their host filament, whereas for halos with M < Mq, the spin tends to be aligned with the closest filament. Right: 
redshift evolution of resp. the filament transition mass (in blue) and the tidal tensor transition mass (in green) derived from the 200 Mpc//i 
ACDM simulations as discussed in the main text. The (cyan) dotted line represents the spin-filament mass transition for a larger smoothing (7.2 
Mpc/h). The displayed error bar is estimated as a 1/3 of the bin mass. The dashed lines correspond to the non-linear masses (for a top-hat filter, 
see Appendix D2) at different sigmas, in particular £7 = 1 (orange) and a = 1.686 (red). The redshift evolution of the transition masses are in 
qualitative agreement with that of Mnl(i7 1.686) though they seem to remain close to power-laws throughout the explored range of redshift. 



scaling, Mq oc Rjj-^. Note that this deper\dence significantly 
depends on redshift. 

The existence and redshift-dependence of the transition 
mass in spin - structure alignment is supported by study- 
ing the halo's spin direction relative to the orientation of the 
large-scale gravitational tidal tensor. The details are given in 
Appendix A, where we find that the more massive halos are 
preferably aligned with the intermediate principal axis of the 
tidal tensor, while smaller halos show a positive alignment 
with the minor axis (which near a filament is the direction in 
which the filament extends). In this approach, our measure- 
ments give for the transition mass: 

M^t « M|)(l + z)-T', 7, = 3 ± 0.3 , M'q~8- 10^^ Mq, (2) 

in good agreement with the skeleton results, equation (1). 
The somewhat larger amplitude of Mq w.r.t Mq can be ex- 
plained by noticing that the tidal tensor associated with the 
gravitational potential smoothed on a scale Rq effectively 
probes larger scales than the smoothed density field itself. 
Fig. 3 shows that if we boost the smoothing scale used to 
define the skeleton to Rq = 7.2 Mpc//i, Ml^■^ and M*^^■^ will 
match very closely (see also Appendix Dl). 

In fig. 3, we also compare Merit (z) to the redshift evolu- 
tion of the mass scale Mnl(z) that corresponds to the fixed 
a{R,z) (defined in Appendix D2). Several values of the vari- 
ance are of interest to track. One is a {R,z} = 1, that formally 
defines the scale of non-linearity. Another is t7-(R,z) = 1.686 
that corresponds to the characteristic mass scale, M*, of col- 
lapsed gravitationally bound halos at redshift z in the spher- 



ical top-hat model. Even though it is clear from fig. 3 that 
both transition masses M^j.j( and Mpj.;j qualitatively match a 
non-linear mass evolution with a 1.686 at low redshifts, at 
high z they still follow a power-law behaviour, while Mnl 
steepens as it probes the steepening power spectrum at ever 
shorter scales. Thus, at high redshifts, the positive alignment 
between the halo's spins and the nearby filaments extends 
to masses that are effectively higher, in terms of the corre- 
sponding characteristic non-linear mass. Although we do not 
have the full quantitative explanation for this effect, it may be 
related to the fact that the filaments at high-z are more pro- 
nounced due to a steeper power spectrum, and are correlated 
with the shear of the surrounding flow more robustly. Note 
that the detection of halos and filaments at these redshifts 
may be a concern for these intermediate resolution simula- 
tions. Whilst we defer a detailed quantitative understanding 
of the redshift evolution of the mass transition, the rest of the 
paper is devoted to explaining the origin of these mass tran- 
sitions. 



3 SPIN INDUCED BY LSS DYNAMICS 

The measurements of section 2.2 strongly suggest that the 
spin direction of dark matter halos is connected to the cosmic 
web. Indeed, these dark halos are embedded in large-scale 
cosmic flows induced by the successive formation of walls, 
filaments and clusters: matter escapes from the voids to the 
walls then to the filaments before flowing along this latter di- 





Figure 5. Trajectory of dark matter particles, colour coded from yellow (early) to Wue (late) via daric red as a function of time. The four panels 
(from top to bottom and left to right) corresponds to different stages of the winding of two walls around a north-east oriented filament. Once 
the dark matter has joined the filament, it heads towards the bottom left part of the panel in the direction of a more massive node of the cosmic 
network. This process is best seen dynamically on http : //www . iap . f r/users/pichon/spin/ 



rection towards the nodes (Zel'Dovich 1970; Bond et al. 1996; 
Pichon et al. 2011). In this framework, let us now argue that 
the first generation of halos are formed during the wind- 
ing of walls around filaments, which provides them with a 
spin parallel to this direction (and whose amplitude is pro- 
portional to the relative impact parameter of the two walls). 
Conversely, the later generations of halos form by mergers 
along the filaments i.e. in the direction parallel to the mean 
flow (as was first pointed out by Sousbie et al. 2008) and 
therefore acquire a spin perpendicular to the filaments (if the 
orbital angular momentum which is converted into spin dur- 
ing the merger dominates). In this scenario, the transition 
in mass measured in fig. 3 in fact reflects a trend in merg- 
ing generation^. Note that this behaviour occurs on several 
scales simultaneously; this multi-scale signal is probed by 



varying the smoothing scale used to define the filaments in 
Appendix Dl. 

As a first check for this hypothesis, the typical distance 
of DM halos from filaments is computed as a function of their 
mass. The more massive halos are typically found closer to 
the filaments (on average at 0.7 Mpc/h) than low-mass halos 
(found on average at 2 Mpc/h from the core of the filaments), 
which is consistent with our assumption because it implies 
that high-mass halos have reached the center of the filaments, 
where they are more likely to merge in the direction of the 
flow. Conversely, low-mass halos (for which large-scale dy- 
namics has been frozen in at an early stage) are found fur- 
ther from the core of the filaments where they are less likely 
to merge. This is qualitatively consistent with fig. 1, which 
shows the distribution of massive halos within a slice of the 
simulation. 



^ Technically it was not possible to build merging trees a posteriori 
on the Horizon 4ti simulation as it would have required too many 
snapshots and thus too much disc space; so we did carry out the 
tests described in Section 3.2 on a smaller simulation to confirm the 
relevance of merger generation as the key parameter. 



3.1 Winding up of dark matter flows 

Let us have a look at fig. 5, which displays the temporal evo- 
lution of dark matter particle trajectories in the vicinity of a 
dark matter filament. This simulation has the special feature 
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that the small-scale modes were erased from the initial con- 
ditions in order to facilitate visualization of the large-scale 
flow, see Pichon et al. (2011) for details. We refer to this sim- 
ulation as a "smoothed" simulation. Here most of the flow is 
in fact embedded in a large wall parallel to the plane of the 
figure. The dots of different colours (from yellow to blue via 
darJc red) correspond to the position of the same dark matter 
particle at successive time steps, and thus allow us to visually 
follow the trajectory of DM particles. The top left panel cor- 
responds to a snapshot somewhat before the flow has signif- 
icantly shell-crossed around the north-east filament. The DM 
particles are plunging towards their filament, while flowing 
within the two walls. On the top right panel, some level of 
shell crossing has occurred in the north east part of the fil- 
ament, and the corresponding particles have started inflect- 
ing their trajectories to wind up around the locus of that fil- 
ament. Since these particles typically have a non zero impact 
parameter relative to the center of mass, as they wind up, 
they convert their orbital motion into spin while generating a 
virtualized structure. Later on, (bottom left panel) this struc- 
ture sinks along the north-east filament towards a more mas- 
sive clump (off field). Meanwhile, the process of DM winding 
from the walls around the main filament continues, and feeds 
(as a yellow trail) the dark matter halo along its current spin 
axis (which is aligned with the axis of the filament). Finally 
on the bottom right panel, another such halo has formed fur- 
ther down the filaments, and we can anticipate that their up- 
coming merger will lead to a structure whose spin's direction 
will be a mixture of their initial spin, and the spin perpendic- 
ular to their relative orbital plane. In section 3.3 below we 
will revisit this scenario using hydrodynamical simulations, 
which allow us to visually identify the spin of forming galax- 
ies. 

3.2 The progenitors of dark halos via merging trees 

In order to understand the previously described mass transi- 
tion and its redshift evolution, we used the code TreeMaker 
(Tweed et al. 2009) to track down the progenitors (both dark 
halos and imresolved flow) of given halos in conjunction 
with their spin orientation relative to the nearest filament. 
TreeMaker involves two steps: first halos are identified using 
a halo finder - in our case a FOF algorithm (Huchra & Geller 
1982; Zeldovich et al. 1982; Davis et al. 1985), while the prop- 
erties of these structures (mass, angular momentum ...) are 
measured. As a second step, the individual DM particles 
which belong to each halo are tracked back in time so as to 
build a merging history tree which regroups all of its progen- 
itors and their properties as a function of time. A relatively 
large number of merging trees were computed and the di- 
rection of the spin of the progenitors relative to their host 
filament was calculated. 

After visual inspection of a subset of those merging 
trees, it was found that: i) the high-mass halos (i.e above 
the critical mass) with a spin perpendicular to their filament 
tend to have a similar history: they often acquire a significant 
amount of mass via a major merger, which is accompanied 
by a significant spin adjustment from a direction initially 
aligned with to a direction mostly perpendicular to their fil- 
ament (and aligned with 62, see Appendix A); ii) in contrast, 
low-mass halos (i.e below the critical mass) are not the result 
of major mergers; often no mergers at all are found at the 



mass resolution of the N-body simulations; those who have 
a spin parallel to the filament seem to have acquired this spin 
direction at a time (the so-called formation time) when they 
have acquired most of their mass by diffuse accretion. 

This behaviour is quite generic as we observed it for a 
few tens of randomly chosen halos. It is illustrated in fig. 6 
where merging trees of two different halo (extracted from 
one of the 50 Mpc/h ACDM simulations) are shown: the 
evolution of the angle between the closest filament and the 
spin of a given halo is plotted as a function of the redshift; 
the colors encode the fraction of mass of the progenitor with 
respect to the final halo mass. The right panel corresponds 
to a low-mass halo (2 ■ lO^^Mij at redshift corresponding 
to more than 300 particles) which forms at redshift z ~ 1.5 
(when it has already acquired more than one half of its mass) 
and suddenly acquires a spin parallel to its closest filament 
at a high redshift (~ 2) close to its formation time. This halo 
does not undergo any significant merger afterwards. Con- 
versely, the left panel provides the merging tree of a high- 
mass halo (8 • W^-^Mq at redshift 0) which forms at lower 
redshift (~ 0.4) as the result of a major merger between two 
less massive halos. This event corresponds exactly to the time 
when it acquires a spin perpendicular to its closest filament. 
What is striking here is the clear flip of the spin direction: the 
two progenitors have a spin aligned with the filament (this 
spin is acquired at higher redshift ~ 1.5) and their merger 
makes the spin of the resulting halo become perpendicular 
to it. These two examples of trees are characteristic of how 
halos below and above the critical mass M^j.jj form and ac- 
quire spin. Note that we actually observe a large dispersion 
of the histories around this mean behaviour. 

To check the statistical robustness of this scenario for 
low-mass objects, let us identify a preferred plane of mo- 
tion at formation time.The following simple test was imple- 
mented: a set of low-mass halos ending with a spin parallel 
to their closest filament are randomly chosen and only ha- 
los for which a time of significant accretion (i.e. their forma- 
tion time) can be determined are retained (this represents at 
least one third of our sample). Their particles are traced back 
one time step before their formation time to quantify the rela- 
tive orientation of their velocities compared to the filament's 
direction. The excess probability distribution of alignment, 
1 + is displayed in fig. 7 and shows that their velocities 
(before formation) are more likely to lie perpendicular to the 
filament (in particular it is foimd that (cos 0y ) ~ 0.47), which 
is in agreement with the scenario (see also fig. Al)). 

In order to assess the statistical relevance of the scenario 
for massive halos (i.e above the critical mass at redshift zero), 
merging trees can also used to determine their most recent 
merger. The mean cosine of the angle between these halos 
and their host filament just after merging is computed and 
compared to the mean cosine of the angle between the closest 
filaments and the spin of the progenitors just before merging. 
This test yields (cos 9) ~ 0.51 before and 0.47 after the last 
merger, which is fully consistent in amplitude with the mean 
cosines found for the Horizon 47r simulation (for which for 
instance (cos0) ~ 0.510 (resp. 0.479) for M ~ IO^^Mg, (resp. 
M ~ IO^Mg)). The 1-a error on the mean is about ±0.02, 
given the size of this sample (only ~ 200 halos massive 
enough in this simulation). This statistical test is in agree- 
ment with the fact that massive halos acquire a spin perpen- 
dicular to their host filament because of mergers. 
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redshift redshift 

Figure 6. Left merging tree of a higli-mass halo (8 ■ W^^M i, at redshift 0). The various colors correspond to different mass fractions (relative to 
2 = 0). The vertical axis correspond to the angle between the host filament and the spin of this halo and horizontal axis to the redshift. Here, 
the spin becomes perpendicular to its filament after an important merger at redshift ^^0.5. Right merging tree of a low-mass halo (2 • 10" Ms at 
redshift 0). Here, the spin becomes suddenly parallel to its host filament when the halo acquires most of its mass by accretion between redshift 
one and two. 
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Figure 7. Excess probability of the velocity-filament alignment just 
before formation time for a sample of 15 low mass halos (M < Merit). 
The average orientation over these 15 halos is plotted with a dot- 
ted line, while the l-a deviation away from that mean is displayed 
in pink. The progenitors' velocities before formation Ues preferen- 
tially in a plane perpendicular to the closest filament (see also figs. 6 
and Al). 



Altogether, we are now able to reconstruct the history 
of spin acquisition by DM halos following the large-scale dy- 
namics. The less massive objects are born by accretion at high 
redshift during the winding of the walls into filaments. This 
process generates a spin aligned with the filaments. Most of 



the halos of low mass at z = are now formed, and their spin 
will not change much because they have already acquired 
most of their mass while future accretion will not be impor- 
tant enough to have a strong impact of their spin direction. 
This behaviour is illustrated on fig. 6 (right panel) where a 
halo forms and acquires a spin parallel to the direction of the 
forming filament, which then remains in the same direction 
during accretion until redshift zero. 

At lower redshift, filaments are collapsing and thus cre- 
ate a flow along their direction in which the more massive 
halos form by major mergers. During this process, these mas- 
sive halos acquire a spin which is the superposition of the 
spin of their progenitors and the orbital spin from the merger. 
As the motion is along the filament, this orbital spin is in 
the plane perpendicular to it. This process was shown on 
fig. 6 (left panel): the spin of the (less massive) progenitors 
are parallel to the filaments and their merger within the fila- 
ments induces a more massive halo whose spin is now per- 
pendicular to the host filament. For these more massive dark 
matter halos, a competition between orbital spin and intrin- 
sic spin during the merger process has just been highlighted 
and explains the resulting orientation. Indeed, their spin are 
not randomly distributed in the plane perpendicular to the 
filaments but are shown to be correlated with one particu- 
lar eigen-direction of the large-scale tidal tensor. This issue is 
fully addressed in Appendix A. 

It is interesting to compare the picture described above 
with Peirani et al. (2004) who focused on the spin magnitude 
(instead of its direction). These authors claimed that spin ac- 
quisition was dominated by merger events rather smooth ac- 
cretion; fig. 6 suggests that this holds true for the spin's di- 
rection of halos more massive than the critical mass M^j.jj as 
well. 
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Figure 8. Distribution of liydrodynamical tracer particles (at z ~ 9, top panels, and z ~ 8 bottom panels). The web-like filamentary structure 
(top left and right panels) of the gas distribution which ends up in the bulge of a unique galaxy at later times is quite intricate, though one main 
wall in which the largest filament is embedded dominates (seen more clearly edge on, in the top right panel). Note the disc-like features with a 
spin parallel to the filament (represented qualitatively as an arrow perpendicular to the disc). Bottom panels: zoom in at later stage, to visualize 
the process of a merger along the filament, before (bottom left panel) and during (bottom right panel) the merger The spin of the merger 
remnant is a combination of the orbital angular momentum and the initial spin of the progenitors; it can therefore depart from the direction of 
the host filament. Note also the ribbon structure of filaments which corresponds to the locus of the second shock. Visual inspection suggests 
these ribbons become broader with time (as predicted by Pichon et al. (2011) as they advect larger and larger amount of angular momenta from 
the outskirt of the gravitational patch) and tend to lie perpendicular to the main wall; as they reach the protogalaxy, they twist rapidly on outer 
shells and build up its outskirt. 



3.3 Visual inspection using hydrodynamics 

For illustrative purposes, let us now turn to a high redshift 
hydrodynamical simulation, as the cold gas that we will ana- 
lyze here follows more closely the caustics of the cosmic web 
than the dark matter, and thus provides a clearer visual im- 
pression of the process of wall winding and spin acquisition 



(in contrast to coUisionless DM, cold gas does not undergo 
shell crossing but shocks and loses degrees of freedom - its 
motion perpendicular to the shock). We are not concerned 
here with how much respective angular momentum gas has 
with respect to dark matter acquired, but only use hydrody- 
namics as a proxy for pinpointing more accurately the loci of 
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Figure 9. Trajectories of tracer particles in the outer region {left panel) and inner region {right panel). On large scales, the gas departs from voids, 
flows along the walls (red arrows) and wind up in filaments along well defined ribbons (green arrows), forming low-mass halos. In the inner 
region of the tracks, the flow is indeed along the filament and spiraling in into a disc whose axis is roughly aligned with the filament. 



shell crossing and identifying the spin axis of galaxies. As we 
will argue later, it is also of interest to consider in parallel the 
environment of low-mass high-redshift and high-mass low- 
redshift halos. Following Dubois et al. (2011) we use tracer 
particles of the gas in a cosmological hydrodynamical simu- 
lation (which is described in Section B) to illustrate this wind- 
ing of walls and the loci and orientations of galaxies. 

Fig. 8 represents the web-like filamentary structure of 
the gas at 2 = 9 in a field of size approximately 50 kpc/fe 
across; the ensemble of tracer particles initially makes up a 
sheet-like structure with a dominant filament embedded in 
it. Note that these tracer particles represent a biased sub- 
set of all tracers as they are chosen so as to end up within 
the bulge of the main galaxy of this zoom simulation at 
some later stage. The tracer particles flow from this sheet 
into the filaments where they form "protogalaxies". The gas 
thereby typically has a non zero impact parameter relative 
to the filament and protogalaxies thus acquire a spin par- 
allel to the filament in which they form (see arrows on the 
figure). These young galaxies then migrate along their fila- 
ments and merge with other galaxies. The spin of the merger 
remnant is a combination of the orbital angular momen- 
tum of the collision and the initial spin of the progenitors; 
it can therefore depart from the direction of the host fila- 
ment. We provide an animation which allows to better see 
this at http : //www . iap . f r/users/pichon/spin/. Fig. Bl il- 
lustrates quantitatively the visual impression of fig. 8 while 
measuring the spin of the circum-galactic discs (between 
O.lRjjir and 0.25Rj„>). 

In fig. 9, a small subset of these tracer particles are ran- 
domly chosen and followed for a while from early (left panel) 
to late times (right panel). On large scales (at early times), we 
note that the flow is indeed dominated by the winding up 
of matter from the main wall around the main filament. The 
tracers' trajectories start perpendicular to the filament within 
the walls. As they reach the filament, they take a sharp turn, 
losing their transverse motion and flowing along the filament 



(left panel). In doing so, since the laminar flows on oppo- 
site sides of the wall will typically have different impact pa- 
rameters, they generate a spinning structure whose axis will 
be aligned with the filament. This spaghetti structure con- 
verges into a quite narrow and elongated plait on either side 
of the forming disc. Given its orientation the induced disc 
will therefore advect secondary infall at its periphery prefer- 
entially along its spin axis (as both the galaxy and its upcom- 
ing secondary infall were assembled via the same winding 
process). 

In a nutshell, the cold gas dynamics of large-scale cos- 
mic flows provides a much clearer illustration of the scenario 
we outlined for the formation of halos with a spin mostly 
aligned with the local filament. We do not discuss any further 
the properties of the alignment of the disc w.r.t. the filamen- 
tary structure given that there exists many caveats (though 
see Section 4.2 below and Appendix B). 



4 CONCLUSIONS & PROSPECTS 

In this paper, the Horizon 47r N-body simulation was used 
to investigate the correlations between the spin of dark mat- 
ter halos and their large-scale environment. For filaments 
defined over a smoothing scale of 5h^^ Mpc, a statistically 
significant signal was detected, indicating that the orienta- 
tion of the spin of dark matter halos is sensitive to the cos- 
mic environment. A mass dependence of this signal was also 
robustly established: low-mass halos are more likely to be 
aligned with large-scale filaments (with an excess probability 
of 15 %) whereas more massive halos tend to be perpendic- 
ular to these (with an excess probability of 12 %). The mass 
transition was found to be redshift dependent and to vary 
like Mcrit(z) « Mo • (1 + z)-T with M'q ~ 5(±1) • IO^^Mq 
and 7s = 2.5 ± 0.2. This critical mass is also found to in- 
crease with the smoothing length (Appendix Dl). These re- 
sults are in agreement with those presented in Appendix A, 
which are derived using the more classic approach of consid- 
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ering excess alignment of the halo spin with the tidal tensor 
eigen-directions. Since the tidal tensor probes larger scales 
than the skeleton of the density, the mass transition is de- 
tected for a larger mass ~ 8(±2) ■ IQI^Mq, which scales 
slightly differently with redshift, i.e. jt = 3 ± 0.3 . Both red- 
shift evolutions are roughly consistent with that of the formal 
non-linear mass scale. 

4.1 Discussion 

This unambiguous result confirms and quantifies recent 
findings (Bailin & Steinmetz 2005; Aragon-Calvo et al. 2007; 
Hahn et al. 2007a; Paz et al. 2008). It is also consistent with 
Sousbie et al. (2009) who found that the spin of dark haloes 
tend to be perpendicular to filaments. Indeed, weighting 
our statistics by the spin magnitude (which corresponds 
to their strategy), we recover their results. Zhang et al. 
(2009) found the same result, but did not have enough 
statistics to probe the signal above the critical mass. The 
study of Faltenbacher et al. (2002) focused on halos above 
lO^'^h^^M-Q, i.e. above the critical mass at redshift zero, 
which is why, in agreement with this work, they found a 
trend for these halos to have a spin perpendicular to their 
closest filament. In contrast, Hatton & Ninin (2001) consid- 
ered halos with mass aroimd 10^^ ^^^M;.,, i.e below the criti- 
cal mass and thus found an alignment of the spins with the 
filaments, again in agreement with this work. This paper also 
confirms the very recent findings of Libeskind et al. (2012) 
who claimed that low-mass haloes (around IO^^^^^Mg:) tend 
to be aligned with the filaments and to lie in the plane of the 
walls. 

Note that Appendix A predicts a strong trend for the 
spin vector of DM halos to lie in the plane of large-scale 
walls. Such a signal was claimed to have been detected in 
observations by Lee & Erdogdu (2007); Trujillo et al. (2006); 
Navarro et al. (2004) and should be re-investigated in view 
of this paper's predictions. Furthermore, observers should 
now also be able to investigate the spin-filament correlations 
in a way directly comparable to the theoretical predictions 
for DM presented in this paper (fig. 3). Indeed, the code Dis- 
PerSE (Sousbie 2011) which identifies filaments using persis- 
tence, can accurately deal with discrete and sparse datasets 
and should provide a good estimator for the direction of ob- 
served filaments. 

The time evolution of the angular momentum of indi- 
vidual dark matter halos obtained by following their pro- 
genitors using merger trees suggests that the spin direc- 
tion results on the one hand from the winding of the walls 
into filaments (first generation, low-mass halos) and on the 
other hand from significant mergers occurring along those 
filaments (second generation, more massive halos). More 
specifically, the arguments we developed throughout the 
paper strongly suggest that the measured correlations can 
be understood as a consequence of the dynamics of large- 
scale cosmic flows. Indeed, low-mass halos mostly form 
at high redshift within the filaments generated by collid- 
ing/collapsing walls. Such a process naturally produces a 
net halo spin parallel to the filaments. In contrast, high-mass 
halos mainly form by merging with other halos along the 
filaments at a later time when the filaments are themselves 
colliding/collapsing. Therefore they acquire a spin which 
is preferentially perpendicular to these filaments. Visual ex- 



amination of 'smoothed' dark matter and hydrodynamical 
simulations lend extra support to this picture (see also Ap- 
pendix B). Measurements of the orientation of the spin rela- 
tive to the eigenvectors of the tidal tensor are also consistent 
with such a scenario provided one takes into account the fact 
that they probe typically larger scales of the density field. 

From the point of view of a large-scale filament, most 
low-mass halos are formed early from patches that are 
part of a planar, flattened inflow of matter onto that fila- 
ment. For gaussian random fields, the tidal tensor in such 
patches is correlated with the filament's direction (via the 
shape parameter 7, Pogosyan et al. 2009), resulting in the 
preferential alignment of the spin of such halos along that 
filament. This process is related to the theoretical predic- 
tions of Pichon & Bernardeau (1999) who demonstrated, us- 
ing the Zel'dovich approximation, how vorticity was gener- 
ated during the first shell-crossing. This vorticity will lie in 
the plane of the forming walls. Extending their predictions 
while focussing now on a 2D flow, we speculate that sec- 
ondary shell-crossing will lead to the formation of vortices 
aligned with the forming filament (see Figures 5 and 7 of 
Pichon & Bernardeau 1999, a possible section perpendicular 
to the axis of the filament)^. In turn these vortices could ac- 
count for the spin of protogalaxies, as was suggested by the 
referee of that paper. 

The excess probability in figs. 3 and A2 lies at the 15-40 % 
level. As such it mainly reflects a residual trend of coherence 
inherited from the large-scale cosmic environment in which 
halos form. This does not preclude the multi-scale hierar- 
chical clustering process to erase part of this more orderly 
dynamics. For instance, clustering and merging on smaller 
scales will in part perturb the large-scale ordered ribbons (as 
is already visible on the bottom panels of fig. 8). 

The trend in fig. 3 is found on the one hand in these mea- 
surements at the high end of the mass function at redshift 
zero (whose dynamics is only mildly non-linear), and on the 
other hand for the gas (a proxy for DM caustics) at very high 
z for lower mass galaxies (again not very far from linear dy- 
namics at this epoch). Hence we found in both regime that 
some imprint of the large-scale structure geometry and dy- 
namics is directly responsible for the spin of the forming ob- 
ject and its post merger transition. In this paper we tried 
and explained the origin of the statistical signal, but we do 
not argue that large-scale dynamics dominates the non-linear 
regime of galaxy evolution on an individual object basis. 

A remaining task involves understanding in details the 
redshift-dependence of the transition masses M^rit described 
by equations (1) and (2) (and also its dependence with the 
smoothing length, which might allow us to identify a scale 
at which the (anti-)alignment is strongest). It would also 
clearly be of interest to further quantify (using larger sam- 
ples) the findings of Hahn et al. (2010) regarding the align- 
ment of the stellar disc/ circum-galactic medium with respect 
to the LSS and investigate how these results depend on sub- 
grid physics and feedback processes. 

Let us conclude this paper by some speculations about 



^ Note however that such a typical caustic should have inherited 
some level of asymmetry (Pichon et al. 2011), which could imply that 
one vortex dominates. 
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what these results imply specifically for the process of galaxy 
formation at high redshift. 

4.2 Implications for galaxy formation 

We are now in a position to speculate about the implica- 
tion of our findings for high redshift galaxy formation. In 
the view of the robust measurements of Sections 2 and Ap- 
pendix A and the visual inspections of Section 3, it appears 
that i) galaxies form preferentially along filaments and ii) that 
their internal dynamics (hence their morphology) inherit im- 
portant features from this anisotropic environment. The first 
point is backed by the distribution of filaments at the virial 
radius (Appendix A of Pichon et al. 2011; Danovich et al. 
2011, and indirectly by fig. 1), the second point by the mea- 
surements reported in this paper. 

Indeed, filaments can be thought of as the loci of con- 
structive interferences from the long-wavelength modes of 
the initial power-spectrum. On top of these modes, construc- 
tive interferences of high frequency modes produce peaks 
which thus get a boost in density that allows them to pass 
the critical threshold necessary to decouple from the over- 
all expansion of the Universe, as envisioned in the spherical 
collapse model (Gunn & Gott 1972). This well known biased 
clustering effect has been invoked to justify the clustering of 
galaxies around the nodes of the cosmic web (White et al. 
1988). It also explains why galaxies form in filaments. In 
walls alone, the actual density boost is not sufficiently large 
to trigger galaxy formation. We therefore argue here that, sta- 
tistically, the main nodes of the cosmic web are where galax- 
ies migrate, not where they form; galaxies generally form 
while reaching filaments from walls. They thus inherit the 
anisotropy of their birth place as spin orientation. During mi- 
gration, they may collide with other galaxies /halos and erase 
part of their birth heritage when converting orbital momen- 
tum into spin via merger. 

Recently, Pichon et al. (2011) showed how the cold gas 
drains out of the prominent voids in the cosmic web, into 
sheets and filaments before it finally gets accreted onto dark 
matter halos. Interestingly, the imprint of the larger-scale 
pancake-structure of the typical cosmic web around a fila- 
ment and a peak allows us to be more specific about the ge- 
ometry of this process. Indeed, one of the striking features 
of figs. 8 and. 9 (probably best seen in the animation online) 
are these ribbon-like caustics which feed the central galaxy 
along its spin axis from both poles. Generically the gas in- 
flow in the frame of the galaxy is double helix-like along 
its spin axis; this is mostly wiped out in DM (and hardly 
visible in fig. 1) because of shell crossing, but quite visible 
for the gas. These ribbons are generated via the same wind- 
ing/folding process as the protogalaxy, and represent the 
dominant source of secondary filamentary infall described 
in Pichon et al. (2011), which feeds the newly formed galaxy 
with gas of well aligned angular momentum (whose direc- 
tion was set by the impact parameter offset of the two neigh- 
bouring walls, which can in turn be attributed to the dissym- 
metry of the four neighbouring voids). As such, the larger- 
scale geometry of the LSS (which biases the formation pro- 
cess) squashes the average neighbourhood of a peak (6 sad- 
dles, 8 voids, 12 walls), into a simpler effective geometry (one 
wall and one embedded filament dominating). Formally, the 
most likely "crystal" of the universe - subject to the con- 



straint of collapse along two axis on larger scales, differs 
from the azimuthally-averaged cubic centered crystal found 
in Pichon et al. (in prep.): it is quite flattened and domi- 
nated by one ridge (Pichon et al. 2011; Danovich et al. 2011; 
Dubois et al. 2011). Note that the gas flowing roughly paral- 
lel to the spin axis of the disc along both directions will typ- 
ically impact the disc's circum-galactic medium and shock 
once more (as it did when it first reached the wall, and then 
the filaments, forming those above mentioned ribbons), radi- 
ating away its vertical momentum (see fig. 10 and Tillson et 
al, in prep.). 

Our speculations here have focussed on a two-scale pro- 
cess. Given the characteristics of ACDM hierarchical clus- 
tering one can anticipate that this process occurs on sev- 
eral nested scales at various epochs - and arguably on var- 
ious scales at the same epoch^. In other words, one expects 
smaller-scale filaments are themselves embedded in larger- 
scale walls (as discussed in Appendix A to reconcile our ex- 
cess alignment with the eigenvectors of the tidal tensor). The 
induced multi-scale anisotropic flow transpires in the scal- 
ing of the transition mass with smoothing, as discussed in 
Appendix Dl. 

Another issue would be to estimate for how long 
this entanglement between the large-scale dynamics and 
the kinematic properties of high redshift pervades. Indeed, 
Ocvirk et al. (2008) have shown that at lower redshift, the 
so-called hot mode of accretion will kick in; how will hot 
flows wash out/ disintegrate these ribbons? Given that they 
locally reflect the large-scale geometry, will the gas continue 
to flow-in along preferred directions (as does the dark matter, 
see e.g. Aubert et al. 2004), or does the hot phase erase any 
anisotropy? Will the above-mentioned smaller-scale non- 
linear dynamics eventually wash out any such trace? 

Finally, note that the actual spin of the stellar disc at 
low redshift need not be trivially related to that of its larger- 
scale gravitational patch (see for instance Hahn et al. 2010), 
as a significant amount of angular momentum redistribution 
takes place in the circum-galactic medium (Kimm et al. 2011) 
over cosmic time. 
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Figure 10. a meridional projection through the velocity flow and gas 
density around the NUT galaxy at redshift 5.1 (Kimmetal. 2011). 
The disc plane is along the horizontal axis, the circle marks the virial 
radius (17.6 kpc). Note the ribbon-like cold-flows seen directly in the 
gas density and the velocity flow, which hit the galactic disc roughly 
along its spin axis. Note also how the gas takes a sharp turn when it 
reaches the ribbons. 

R Bouchet for allowing us to use the magiqueS supercom- 
puter during commissioning, and to S. Rouberol for making 
it possible. We also thank D. Mimro for freely distributing his 
Yorick programming language and opengl interface (avail- 
able at http : //yorick. sourcef orge .net/). 
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APPENDIX A: SPIN-TIDAL TENSOR CORRELATIONS 

Let us present a complementary set of measurements: 
the correlations between the spin axis of dark halos, and 
the orientation of the large-scale gravitational tidal tensor 
Tjj = djjcp — ^AipSij. To describe the orientation of the tidal 
tensor, we define ei, 62 and 63 to be the minor, intermedi- 
ate and major eigen-directions of T,y according to the sorted 
eigenvalues Aj ^ A2 ^ A3 of the Hessian of the gravita- 
tional potential, (with which the tidal tensor shares the 
eigendirections) . 

Besides works based on the correlations between spin 
orientation and the cosmic web described in the main text 
(Section 1), the only numerical study to date looking at the 
alignment between halo spin and tidal tensor was done by 
Porciani et al. (2002) who predicted its orthogonality with the 
major principal axis but also found that galactic spins must 
have lost their initial alignment with the tidal tensor pre- 
dicted by TTT. Direct observations of the alignment between 
the spin and the tidal tensor eigenvectors have also been car- 
ried out: the first attempt by Lee & Pen (2002) studied the 
correlations between the disk orientation of the galaxies from 
the TuUy catalog and the shear reconstructed from The Point 
Source Catalog Redshift survey and rejected the hypothe- 
sis of randomness at a 99.98 % confidence level. More re- 
cently, Lee & Erdogdu (2007) detected some correlations be- 
tween the spin and the intermediate eigenvector of the tidal 
tensor and found that galactic spins were also preferentially 
perpendicular to the major principal axis but this signal re- 
mains weak. To overcome this lack of a clear numerical de- 
tection, the use of the 43 million halo sample of the Horizon 
47r simulation presents a tremendous advantage, as it allows 
us to very robustly calculate the correlations between halo 
spin orientations and the local tidal tensor. We show in this 




Figure Al. The walls of the potential (in blue) and its filaments (in 
cyan), together with the velocities of dark matter particles (in red) 
of a 20 Mpc/h 512^ dark matter particle ACDM simulation with 
WMAP-1 cosmogony. At redshift zero, most dark matter halos sit 
in these walls, while the velocity field empties the voids and flows 
within those walls. This divergent flow is best seen in the top (hence 
bottom) left void (see also the animations online). 



Appendix that these measurements are not only consistent 
with the spin-filament correlations, but actually lend addi- 
tional support to the interpretation of this paper in terms of 
large-scale dynamics along walls and filaments. 

At the onset of non-linearity, the gravitational potential 
tracks the velocity potential of the matter flow. Thus, the 
signs of A,'s determine whether the flow in the correspond- 
ing direction compresses (A; > 0) or rarities (A, < 0) matter. 
At a given smoothing scale, this criterion can be used to par- 
tition space into peak-like < Aj ^ A2 ^ A3, filament-like 
Al < 0, < A2 ^ A3, wall-like Aj < A2 < 0, < A3 or void- 
like Aj ^ A2 ^ A3 ^ 0, regions (Pogosyan et al. 1998). From 
this point of view, in the peak regions matter compression 
is strongest along 63 and weakest along Cj. In the filamen- 
tary regions, Cj gives the direction of the filament, while the 
walls are collapsing along 63 and extend, locally, in the plane, 
spanned by Cj and 62. 

At this stage, it is important to note that the tidal field 
probes larger-scale structures than the filamentary struc- 
ture studied in Section 3 as the gravitational potential is 
a smoother version of the density field (through Poisson's 
equation). In other words, the skeleton of the potential 
(which locally corresponds to the eigen-directions of the tidal 
tensor) traces the cosmic structures (walls, filaments,. . . ) on 
scales much larger than the skeleton of the density. Thus, if 
we turn to a formulation in terms of the skeleton of the po- 
tential, the filaments described in Section 3 are embedded 
in the large-scale walls of the potential field (as illustrated 
in fig. Al), therefore protogalaxy formation begins with the 
first collapse (namely the collapse of 63) leading to the forma- 
tion of the large-scale walls because it corresponds on smaller 
scale to the time when the filaments (of the density field) 
form by winding. 

Note that a stricter definition of a filament in the local 
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theory of the skeleton (Sousbie et al. 2008) as a ridge in the 
density profile, associates its local direction with the minor 
eigen-direction of the Hessian of the density, (here, to match 
enumeration, taken with a negative sign) . Thus, as the po- 
tential is two derivatives away from the density, the excess 
probability of alignment between the halo spin and Cj should 
be quite similar to that of the alignment between the halo 
spin and the filament's direction. 



Al Alignment between spin and tidal eigen-directions 

In order to compute the excess probability of alignment be- 
tween the spin of dark matter halos and the eigenvectors 
of the tidal tensor, the density field is again smoothed over 
5 ?!^^Mpc (for z=0), Poisson's equation is solved via Fast 
Fourier Transform in smaller overlapping boxes (boundary 
effects are foimd to be insignificant inside the boxes), the 
Hessian matrix of the potential is computed using a finite 
difference scheme ans this matrix is finally interpolated to 
halo positions. We then measure the angle between the an- 
gular momentum vector of the halo and each eigen-direction 
of the tidal tensor and compute the histogram of the absolute 
value of the cosine of this angle; after normalization, it gives 
1 + f , the excess probability of alignment between the spin 
and the tidal tensor eigen-directions. As in Section 2.2, the 
data is then split by halo mass. 

Fig. A2 displays these excess probabilities in three pan- 
els corresponding to the orientation of the spin w.r.t ei, 62 
and 63 (top, middle and bottom panel, respectively), at z = 0. 
Halos of all masses have spins that preferentially avoid the 
direction of the strongest large-scale compression, 63, at 10% 
excess probability (boffom panel). The exception is the very 
highest mass bin, which shows additional halo population 
with spins aligned with 63. 

For the high-mass, M > M(, w 8(±2) ■ lO^^M©, halos 
we detect a strong trend for the spin to be aligned with 62, 
the intermediate principal axis of the tidal tensor with an ex- 
cess probability of up to 40% (bJue and green lines, middle 
panel), and to be perpendicular to the minor, ei, principal 
axis with an excess probability of up to 20% (bJue and green 
Jines, top panel), a result in agreement with Lee & Erdogdu 
(2007); Porciani et al. (2002). 

Spins of the lower-mass, M < Mq, halos tend, in con- 
trast, to align with ei (red and orange lines, top panel) thus, 
preferring the direction of the filamentary structures ei, with 
an excess probability of up to 15% and in a weaker way to 
align with 62 (with an excess probability of 5%, see red and 
orange lines in the middle panel). 

These results are in exact agreement with our findings 
using the skeleton that the spins of sufficiently large halos 
prefer to be perpendicular to the filament's direction while 
small halos show a positive correlation for a spin orientation 
along the filaments. Indeed, the bottom panel of fig. A2 is 
almost identical to the spin-filament correlation foimd in fig. 
3. For the tidal tensor, the transition occurs at a somewhat 
higher mass, Mq w 8 • lO^^Mr.:, than for the skeleton probe. 
This is not surprising, because of the effectively larger scales 
probed by the tidal tensor and the observation (see Section 
2.2 and Appendix Dl) that the critical mass increases with 
the smoothing length used to define the large-scale structure. 

The redshift-dependence of the transition mass was also 
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Figure A2. Excess probability of alignment between the spin and 
the minor/intermediate/ major axis {from top to bottom) of the tidal 
tensor in the Horizon 4n simulation. Different mass bins are color- 
coded from lO^^(red) to 10^* Mq (blue). A transition is detected: spin 
of high-mass halos tend to be aligned with the intermediate principal 
axis {middle panel), whereas low-mass halos spins are more likely to 
point along the minor axis (top panel). 
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investigated (as described in the main text, see Section 2.3) 
and is found to be: 

M^t w M'q{1 + z)-'^>, 7, = 3 ± 0.3 , M'q- 8(±2) • lO^^Mg. 

A2 Consistency with LSS cosmic flows 

Let us describe this spin acquisition in the framework of 
large-scale dynamics and the ellipsoidal collapse model 
(Lynden-Bell 1964; Zel'Dovich 1970; Icke 1973; White & Silk 
1979; Peebles 1980; Lemson 1993; Bond & Myers 1996; 
Sheth et al. 2001; Desjacques 2008, among others). Low-mass 
halos form by accretion at high redshift. At this time, sheets 
and filaments are forming by the successive collapses of 
(resp.) 83 and 62- When 63 is collapsing, low-mass halos form 
by accretir\g particles whose motion is along the direction 83. 
This process induces a spin perpendicular to this directior\, 
i.e. in the plane (81, 82) which is the plane of the large-scale 
wall in which they are located. Then, 82 begins to collapse 
and other low-mass halos form from objects moving in the 
plane (82, 83 ) acquiring therefore a spin perpendicular to this 
plane i.e. aligned with 81 (which is the direction of the form- 
ing large-scale filament). The first generation of dark matter 
halos (of typically low mass) is now formed; as mentioned 
in Section 3, for low mass halos at redshift zero, the spin will 
not change much as they have already acquired most of their 
mass. Fig. 6 (left panel) provides a clear illustration for this: 
a low-mass halo forms and acquires a spin aligned with 81 
as expected. It has been pointed out above that halo spin 
orientation must be either aligned with 81 (correlated to the 
filament's direction), or in the plane (81,82) depending on 
the time at which they form. This is in good agreement with 
the top and middle panels of fig. A2 which shows an excess 
probability for their spin to be aligned with 81 (or with the fil- 
aments) and in a weaker way with 82 (see the red and orange 
lines which represent the low-mass halos at redshift zero i.e. 
the halos with a mass between 3 • lO" and 3 • W^^Mc,). 

Later (as described in Section 3 in terms of flows along 
the filaments), 81 collapses and the halos located in the fil- 
aments of the potential stream along this direction. Most of 
the more massive halos then form by mergers in this flow 
and therefore acquire a spin which combines the spin of their 
progenitors and the orbital spin provided by the merger. The 
orbital spin must be in the plane perpendicular to 8^ i.e. 
(82, 83) because the progenitors move along the filaments be- 
fore merging, whereas the spin of their (less massive) pro- 
genitors is in the plane (81,82). The resulting angular mo- 
mentum is therefore a superposition of these various spins, 
statistically more likely to be aligned with 82, which is what 
was shown on fig. 6 (right panel) and which is also in good 
agreement with the middle panel of fig. A2, where the blue 
and green lines representing halos above 2 ■ 10^^ Mq re- 
veal a strong trend for these high-mass halos to be aligned 
with 82. For these massive dark matter halos, the competi- 
tion between the orbital spin and the intrinsic spin during 
the merger process was already pointed out in Section 3. The 
excess probability for their spin to be aligned with 82 sug- 
gests that neither one nor the other dominates. Nevertheless, 
very massive halos (with masses above 10^*Mf ,) represented 
with a blue line in Fig A2 seem to have their spin less perpen- 
dicular to 83: actually we can observe two modes, one per- 
pendicular and one aligned with this direction, which can be 



understood if they are the result of a further generation of 
mergers whose intrinsic spins were already perpendicular to 
the filament. 

Let us emphasize that this explanation and that which 
was presented in Section 3 are consistent. Indeed, in Sec- 
tion 3 our claim is that winding of the walls is responsi- 
ble for the direction of the spin of low-mass halos. Mean- 
while, in this Appendix, the focus is on the first collapse 
along 83 i.e. on the formation of walls. However, as we was 
pointed out in the introduction of this Appendix, the tidal 
field probes larger scale structures than the filaments of the 
density field studied in Section 3. Taking into account this 
difference, it turns out that the two analyses are comple- 
mentary (describing the large-scale dynamics on different 
scales). It also helps reconciling the findings of Porciani et al. 
(2002); Lee & Erdogdu (2007) which rely on the tidal tensor 
with those of Bailin & Steinmetz (2005); Aragon-Calvo et al. 
(2007); Hahn et al. (2007a); Paz et al. (2008) which involve the 
density field filamentary structure. 



APPENDIX B: CIRCUM-GALACTIC MEDIUM SPIN 

Let us first describe briefly the hydrodynamical simulation 
that we will use to assess the alignment of the gas compo- 
nent surrounding the galaxies with the large-scale filamen- 
tary structure. This simulation is described in more details 
in Dubois et al. (2011) and corresponds to their SHhr run. 
Let us recall here its basic properties. The SHhr simulation 
follows the formation of a massive (Mvir = 5 ■ 10" Mq at 
2 = 6) halo within a re-simulated region in a 100 Mpc//i box 
size, its dark matter mass resolution is 1.3 ■ 10^ M©, and min- 
imum cell size is 17 pc. The cosmology employed in that run 
is slightly different from the parameters of the Horizon in 
simulation, and is compatible with the WMAP7 cosmology 
(Komatsu2011): Om = 0.27, Oa = 0.73, Of, = 0.045, Hq = 70 
km s^^ Mpc^^, n = 0.961, and og = 0.8. The gas is allowed 
to cool down radiatively down to Tg = 100 K, assuming 
an initial metal enrichment of lO^'^ Z 5 (Sutherland & Dopita 
1993). A UV background heating source term is added to 
the gas energy equation following Haardt & Madau (1996) 
with reionization taking place at Zj-gion = 8.5. Star forma- 
tion is allowed in gas density regions above wq = 50 H • 
cm^'' using a Poisson random process (Rasera & Teyssier 
2006; Dubois & Teyssier 2008) that reproduces the Schmidt- 
Kennicutt law = e^p/fff, where is the star formation 
rate density, = 0.01 the star formation efficiency, and tg 
the local free fall time of the gas with local density p. No 
feedback from supernovae or Active Galactic Nuclei is ac- 
counted for. Tracer particles that passively follow the motion 
of the gas are scattered in the initial conditions and allow us 
to trace back the Lagrangian trajectories of gas elements that 
end up in collapsed structures. 

Fig. Bl displays two large-scale views of this hydro- 
dynamical simulation and its tracer particles at redshift 9 
shown in fig. 8. Structures and substructures in dark mat- 
ter are detected with the Most-massive Sub-node Method 
(Tweed et al. 2009) and only masses above 5 • IO^Mq are se- 
lected. The spin of the circum-galactic medium (accounting 
for non-star forming gas only with gas density below «h < 
50 H.cm^^) between 0.1Rp,> and 0.25Ryir is then computed 
and its orientation is represented with dark red segments. 
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Figure Bl. Distribution of hydrodynamical tracer particles (blue) at 
redshift 9, along two different viewing angles. The spin direction of 
tfie circum-galactic medium between O.lRj,,,. and 0.25R;,/,. of some 
protogalaxies is plotted in red. On average the spins are more likely 
to be aligned with the filaments (as seen with protogalaxies 2,3,4 and 
5 for example). 



This somewhat ad-hoc criterion used to define the spin re- 
flects our focus on the angular momenhim of the secondary 
infall gas, which has just been or is being accreted; our mea- 
surements correspond to an average of the spin within this 
subregion of the dark matter halo. This figure shows a good 
alignment of the spins with the circum-galactic polar fila- 
ments (i.e. the filament which visually flows along the polar 
axis of the galaxy, in particular for the clumps 2,3,4,5). This 
is consistent with visual inspection (fig. 8, but best carried 
online), and with the prediction of fig. 3, since the critical 
mass at that redshift is 6 ■ lO^'^M;:,. The spin of a couple of 
low-mass clumps in that field (noticeably clump 1 and 6), 
are in fact poorly estimated automatically, as the dark matter 
clump center can be offset at that redshift relative to that of 
the circum-galactic disc. 



APPENDIX C: ROBUSTNESS OF SPIN-FILAMENT 
CORRELATION 

Let us assess how robust the spin orientation - filament cor- 
relation found in Section 2.2 is by carrying a few consistency 
tests. 

In order to check the effect of mass resolution (the spin of 
very low-mass halos is poorly defined for instance, as too few 
particles are involved in its measurement), the same mea- 
surement are carried out in a smaller simulation (256'' par- 
ticles in a 50 /i^^Mpc periodic box with the same cosmology) 
for which same physical masses are represented by higher 
numbers of particles: mass bins from 3 • 10^^ (red) to 6 ■ 10^^ 
Mij: (green) correspond to 30-1400 particles in Horizon in 
(corresponding to a lower threshold for the FoF detection) 
since the mass per particle is 7 ■ 10^ M(.,; in contrast, 350 to 
15000 particles are found in halos of the same mass for the 
small simulation since the mass per particle is 6 • 10^ Mq. 
The detection of the the same phase transition (see fig. CI) 
occurring at the same halo mass demonstrates that the signal 
is not induced by limited mass resolution. This result is also 
consistent with the findings of Aragon-Calvo et al. (2007); 
Hahn et al. (2007a). 

Another simple check involves varying the procedure 
by choosing for each halo the closest segments of the skeleton 
(instead of for each segment, the closest halos). The signal we 
get is very similar to fig. 3 (low-mass halos tend to be parallel 
to the filaments with an excess probability of 15% and high- 
mass halos perpendicular with an excess probability reach- 
ing 20% for the more massive bin which is even stronger than 
in fig. 3), suggesting that the measured correlations are inde- 
pendent from the detailed procedure implemented to iden- 
tify neighbours. 

Halos at the nodes of the skeleton cannot have a well 
defined closest-segment direction, as more than one skele- 
ton segment typically qualifies, and could therefore bias our 
measurements. So as to quantify this effect, the same algo- 
rithm is implemented but with a new criterion: halos closer 
than a certain distance to the nodes are not considered. An 
even stronger signal is detected, which leads us to conclude 
that nodes introduce extra noise and are not the cause for the 
observed signal. 

One might also think this result could depend on a den- 
sity threshold for the underlying filament. Hence, the same 
data in low, middle and high-density filaments were plotted: 
the excess probability of spin-filament alignment is found to 
be the same whatever the density inside the filaments is. All 
these tests demonstrate the overall robustness of the mass- 
dependent transition of the relative alignment. 

APPENDIX D: CHARACTERISTIC MASSES 

In the main text the density field was smoothed over a scale 
of 5 Mpc//i corresponding to a mass of 1.9 x lO^^M,:). The 
transition mass foimd in Sec. 2.2 is therefore defined rela- 
tively to this mass. In the context of hierarchical clustering, 
as long as the field is smoothed on scales where filaments are 
still well defined, one can anticipate some scaling of this tran- 
sition mass with smoothing. This transition mass should re- 
flect the connection between the geometry of the larger scale 
flow and the mass scale corresponding to galaxies forming 
and drifting on this cosmic web. 
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Figure CI. Excess probability of spin-filament alignment in various 
simulations. Measurements in the Horizon 47r simulation are plot- 
ted with dotted lines and measurements in the smaller simulations 
are plotted with solid lines. Different colors correspond to different 
halo mass bins from 3 • lO^^(red) to 6 ■ 10^^ M t, (green). The signal in 
the Horizon An simulation is statistically consistent with that of the 
smaller simulations. 



2x10' 



1 



1x10' 



5x10' 



* • • 

Tidal Tensor 
z=0 

Filament 
z=0 



Filament 
z=2 ,»' 



1.5 



2.0 



3.0 

-/^smooth [Mpc/h] 



5.0 



7.0 



Figure Dl. Evolution of the critical mass as a function of the smooth- 
ing length for the filaments at redshift (cyan) and 2 (grey) and the 
tidal tensor at redshift (purple). The critical mass increases with the 
smoothing length with a redshift-dependent slope. The blue dot was 
measured in the in simulation for the filaments. 



Dl Smoothing-dependence of the critical mass 

Fig. Dl displays the evolution of the critical mass with the 
smoothing length. It is found that M^j.;((R) and Mjj.jf (R) can 
be well fitted by power-laws, namely: 



with P{k) the powerspectrum, and the top-hat filter defined 



Mo{R) ~ Mo(-Ro) 



(Dl) 



where Mg(Ro) - 3.6(±1) ■ IQi^, oc' ~ 0.8 ± 0.1, M^(Ro) - 
4.8(±1) • lO^^ and a' ^ 0.28 ± 0.04. This scaling is not in- 
consistent with the discussion of the origin of the alignment 
given in Section 3 in as much as a smoothing length de- 
fines a set of filaments and therefore picks out a halo mass 
scale corresponding to the halos which are flowing/ merging 
along these filaments. Note that in practice the dependence 
on smoothing is actually rather weak (expressed in terms of 
mass, we have M^j-it ^smooth filaments defined via the 
skeleton and Merit ^smooth structure anisotropy is de- 
scribed via the shear tensor). Fig. Dl also shows that a de- 
pends on redshift. Note that as expected, the critical mass for 
the tidal tensor matches that of the filaments smoothed on 
a larger scale (~ 7.5 Mpc/h instead of 5). This confirms the 
idea that the potential is close to a smoother version of the 
density field. 

D2 Non-linear mass evolution 



It is of interest to compare the transition masses, M^j.jj(z) and 
M^j.jj(z) to the mass scale that tracks the development of non- 
linearity in structure formation. The variance of the density 
field smoothed on scale R obeys 



o-^{R,z) = D{zY / P{k)W^{kR)d?k, 



(D2) 



byW2(x) =9(sinx/x- 
is given by 

D(z) = ln,„HlH{z) 



■cosx) /x . The growth factor, D(z) 



(l+z)dz 



H(z)3 



with H(z) = Ho V^X^T+zp + nX- Here Hq, n,„ and Qa 
are the Hubble constant, the dark matter and the dark energy 
density parameters at z = respectively. 

Fixing the level of (non-)linearity by the condition 
a{R{z),z) = const implicitly defines the redshift evolution 
of the smoothing scale R{z) (expressed in comoving Mpc) 
that maintains this level of (non-)linearity. This, in turn, cor- 
responds to the mass scale 



Mnl(z) = ^npR{zf, 



(D3) 



where p is the present-day average density of matter in the 
Universe. 

For the matter dominated CDM Universe with a scale- 
free power spectrum P{k) oc k'\ Mp,jl(z) a (1 +z)^^/'"+^). 
In the Universe with realistic parameters, the redshift depen- 
dence of Misjl(z) is not a power law, both due to the influ- 
ence of the A-term that slows the growth of the structure at 
low redshifts, and steepening of the spectrum as one moves 
to smaller scale at high redshifts. 



